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We present a model for radiative energy loss of heavy quarks in quark gluon plasma which incorporates coherence 
effects. We then study its consequences on the radiation spectra as well as on the nuclear modification factor of open 
heavy mesons produced in ultrarelativistic heavy ion collisions. 
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1. Introduction 

Ultimately, one expects to probe fundamental properties of the hot medium formed in ultrarelativistic heavy ions 
collisions (URHIC) thanks to the jet quenching observed in these URHIC. In this respect, the quenching of heavy 
quarks (HQ) enriches the analysis and permit to test some general concepts like mass-hierarchy 1 1 1. Understanding 
the nature of HQ energy loss is a rather broad and complex subject reviewed f.i. in O. In O, we have advocated that 
the gluon formation length characterizing the radiative energy loss is small enough at RHIC energies to evaluate this 
quantity through a generalization of the Gunion-Bertsch (GB) approach |4| for finite mass and have obtained good 
agreement with the data. As this ansatz is questionable at LHC, we have developed a model to study the impact of 
coherence on HQ quenching and present here our main results. We insist that all shocks at the origin of the radiation 
are treated in a running coupling approach |5 1. Therefore, the large value of as at small momentum transfer prevents, 
in our view, the direct application of most of the available frames designed in recent years for dealing with this topic. 

2. Estimates for the formation times of radiated gluons 

As the regime pertinent to the gluon radiation depends on the comparison between scales pertaining to the medium 
and the gluon formation length /y^, we start by estimating this quantity. In the case of a gluon emitted through a single 
scattering. If can be estimated from the virtuality of various Feynman diagrams contributing to the radiation \i4i . 
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Figure 1 : 3-gluons diagram of the x + g + Q' +g process (left) and typical diagram describing gluon radiation induced in a multiple scattering 
process (right). 
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For the case of the "3-gluons" diagram on fig. [T](left), specific of QCD and contributing to the radiation at mid- 
rapidity, one finds tf = 2 where ^pre = ^ - ^ is the 4-momentum of the pre-gluon and nig is the thermal gluon 

I "-pre 

mass |6i used in (TJ for the evaluation of radiative energy loss along a finite opacity expansion. For 7 » 1 , one obtains 

I'/' ^tf. ^-^^^ ^, (1) 

where M is the quark mass and x the energy fraction of the gluon x = co/E. When M = rrig = q^,we recover the 
usual law // ~ |r 19J. As most of the shocks happen with ^ [i (the Debye screening-mass associated to the 

2{l-x)xE 



medium) and lead to the radiation of a gluon with ~ as well, one gets (//,sing) ~ x^M^+{i-x)ml+^^ ' ^^^^ rUg - jd for 
practical purpose. At x Xcr = ^, // admits a maximum max(//) ^ Accordingly, coherence will aff'ect mostly 
gluons radiated at x ~ and its impact will be reduced for HQ due to the hierarchy //(M^) > //(M^) > //(M/^). 
Comparing If with the mean free path X, one obtains a first qualitative criterion that permits to estimate whether 
coherence eff'ect should be included. For T = 0.25 GeV and as = 0.3, one finds that it is the case for Ec > 10 GeV 
and for Et > 25 GeV. 

We now turn to the estimation of // in a multi- scattering process happening inside a medium considered as infinite. 
Following the semi-classical argument of LP |8 |, we proceed by considering the phase 0(0 = ^ P{t') • k(f)^ 
associated to the radiation process, where P and k represent the instantaneous momentum of the parton and of the 
radiated gluon. Such a phase is also found in diagrammatic approach stemming from QFT |9 |, up to recoil corrections 
for finite x. Assuming that the rescattering primarily aff'ects the gluon due to its larger color charge, one finds 



<o(/)) ^ I 



£2 co^ co^2 



I p 

(2) 



»^/,sing t ^ 



/,mult 



for E M, where q is the gluon transport coeflftcient, that is the mean square momentum transfer per mean free patlj^ 
Ag. The second term in the r.h.s reflects the contribution of the rescatterings to the increase of the phase, that is to the 
separation of the gluon from its radiating parton. // is then defined as the length for which O = Odec ~ 1- Accordingly, 
one gets // ^ niin(/y^^sing, //,muit — 2 ^^/O^^^ZV§). First, we note that the coherence eff'ects mainly show up when the 
second term dominates over the first, that is for > L** = ^ — . For finite x and M nig, L** » A, 

J ^ vising q * 

meaning that coherence eff'ects do not aff'ect significantly the spectrum for // g [/1,L**]. Comparing If^^mg and //,muit, 
one obtains the following regime^ already identified in O : a) (low energy regime) for E < E^q lpm ~ Mm^ Iq, 
^/,sing < //,muit and coherence eff'ects are thus negligible; b) (intermediate energy regime) when E > £^noLPM, one 
has //,muit ^ ^/,sing for X e [xi, X2] with xi ~ mpqE and X2 ~ -^qE^/M^ and coherence eff'ects should be taken into 
account on this domain; c) (high energy regime) when E > £'lpm ~ M^/q, X2 ~ 1 and coherence eff'ects dominate 
the full spectrum, except for a small domain at 61; < A{^^±} ~ 7". In this large energy regime, the formation length 
becomes mass-independent and one expects similar radiation for all flavors. 



3. Model for gluon radiation induced by multiple collisions 



1 ^.(icrgQ^gQg 



Let us recall that one obtains the power spectrum of gluons induced by single scatterings as the ratio ^ 
that is following relation (12) of 131: 



d(T 



dNg 2N,a,(l - x) C . (, , gj U , , . , ^hT ... 



el 



is the probability density associated with a single (Rutherford) scattering. If we now consider the gluon radiated 
after some rescattering by the medium, as illustrated in fig. [T] (right), our basic hypothesis is that the quark - gluon 

^Including diagrams where the quark undergoes rescattering leads to an effective MFP X = 2Ag which should be taken to calculate the q. 
^Further discussion can be found in flOl . 



Hard Probes PB Gossiaux / Nuclear Physics A 00 (2012) 7-g| 



3 



system cannot resolve the various kicks qi, q2, . . .qn during the formation length, so that all scattering centers act as 
an effective single one shocking the system with a transverse momentum = + • + q_Ln- Hence, we could extend 
eq. d3b by substituting Pi(q±) ^u^n^oSQa.)^ where A/^coh = j the average "coherence number" and !P„ is obtained 
frornPi through convolution. However, it is important to note that the multi-scattering process causes the radiation to 
happen in a shorter formation length and cures the collinear divergency that would appear for massless partons. Those 
two facts are correlated as the amplitude of the gluon field is oc If [11 J . so that the diff'erential power spectrum dl oc ij,. 
More precisely, the finite-mass GB transition probability writes 



iMoBr oc 



I /,sing ^ 



- 4±)\ 



(4) 



so that the denominator in the In of eq. ^ can be seen as a reminiscent of the collinear formation length in our finite 
gluon mass approach: -h x^M^ = jp^. Based on this interpretation, we propose to replace the denominator in the 

In by a quantity evolving smoothly from -jj^ > when coherence eff'ects set in: 



m i -h X M = 
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m„ = 
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= ml + x^M^ + ^{qi) = ml + x^M^ + 



2<D, 



dec 



qoj 



(5) 



Performing these changes in eq. ^ and dividing by the formation length // = NcohX one obtains an eff'ective power 
spectrum per unit length at small x, valid for small or large Ncoh- 



j2 7X^1 

^ ^ model 

dzdoj 



(6) 



Starting from this expression, one can recover all regimes identified in section |2] In particular, for Ncoh > 1 but not 

too large and <^^) ^ m^, one finds J-^ ^In^l -h = 0(N^^^) which indicates that coherence eff'ects do not alter 

the spectrum significantly. This corresponds to the case A < If < L^^ discussed above. To calibrate our model, we 
consider the analytical results established in |[T2ll for = M = and A^coh » 1 • 
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where k = 



2(jj 



(J) 



To make the connexion to this result, one has to consider more seriously the case of Coulomb scatterings, for which 
^1 ^ ( 2^ 2^2 - For a given A^coh » 1, one finds O that {Q\) oc N^oh V^cohy^^ instead of (2^) oc NcohJ^^ for Gaussian 

diff'usion. This leads to an eff'ective q = ^ x ln^2,^^^j and a reduction of the formation length If 1^^^ = 
/// ^In ^2 -yj^^^^^ j in our expression (6 ) which then shows the ^^In I dependence of eq. (7 ), with 
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(8) 



dec 



The first factor reflects the diff'erence between the color factor associated to the planar diagrams that dominates in the 
large A^^ large A^coh limit 1 12] and the color factors in a single process while the second reflects the "quality" of our 
efl'ective scattering center approach. For practical purposes, we choose Odec = 2 and rescale our model (l6|) by a factor 
= -jf- when A^coh » 1, in order to guarantee quantitative agreement with d7b; we also evaluate Pi as well as 
the transport coefficient q entering our expressions with the one eff'ective gluon exchange model presented in f5l . In 
fig. [2] we display some radiation spectra per unit length resulting from our model for typical parameters encountered 
in URHIC. 
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Figure 2: Power spectra per unit length for gluon emission through a single scattering (GB) and for gluon emission through a multiple scattering 
process modeled by a single effective scatterer (LPM); see text for details. 



For /?-quarks, 20 GeV/c ^ ^^noLPM- Thus, one is still in the low energy regime and impact of coherence on the 
spectrum is pretty weak; for c-quarks, 10 GeV/c and 20 GeV/c belong to the intermediate energy regime. The effect 
is more pronounced and increases with E as expected. 



4. Consequences for heavy quark quenching in ultrarelativistic heavy ion colhsions and conclusion 

In fig. [3] we present the nuclear modification factor (Raa) of Z)-mesons in central Au-Au (RHIC) and Pb-Pb (LHC) 
collisions. The "El. + rad. GB" curves correspond to the calculation presented in |3 |, with as = 0.3 in eq. Q and a 
rescaling of the interaction rate by a factor K = 0.7 introduced to reproduce the Raa of non photonic single electrons 
measured at RHIC. The "El. + rad. LPM" curves correspond to the model presented in section [3] One observes a 
typical increase of 40% (resp. 100%) due to coherence for the Raa(D) at pr = 10 GeV/c (resp. pr = 50 GeV/c). 
This originates from the reduced quenching in the "rad. LPM" case (see fig. |2]), which is also responsible for the 
increase of Raa with pr for ^ 15 GeV/c. We conclude that the pr-vsinge available at LHC off'ers better chance to 
discriminate between models of HQ energy loss, although it will require good experimental accuracy. 
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Figure 3: Dashed: Raa of £)-mesons obtained with our MCof^HQ generator 1 5 1 with elastic + radiative energy loss evaluated by extending the GB 
approach to massive quarks; plain: same with coherence effect included for radiative energy loss. 
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